turn and classical aspects, and therefore, 
it will be useful for implementation of 
quantum computing. 

Recent advances in quantum informa- 
tion theory have inspired an explosion 
of interest in new quantum algorithms 
for solving hard computational prob- 
lems. Three basic “non-classical” proper- 
ties of quantum mechanics — superposi- 
tion, entanglement, and direct 
tensor-product decomposability — were 
main reasons for optimism about capa- 
bilities of quantum computers and 


quantum communications as well as for 
a new approach to cryptography. How- 
ever, one major problem is keeping the 
components of a quantum computer in 
a coherent state, as the slightest interac- 
tion with the external world would cause 
the system to decohere. Another prob- 
lem is measurement: by the laws of quan- 
tum mechanics, a measurement yields a 
random and incomplete answer, and it 
destroys the stored state. 

This proposed reinterpretation of 
quantum formalism opens up new ad- 


vantages of quantum computers: if the 
Madelung equations are implemented 
on a classical scale (using, for instance, 
electrical circuits or optical devices) , all 
the quantum effects important for com- 
putations would be preserved; at the 
same time, the problems associated with 
decoherence and measurement would 
be removed. 

This work was done by Michail Zak of Caltech 
for NASA’s Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 
NPO-46731 


© Optimal Padding for the Two-Dimensional Fast Fourier Transform 

Appending data to an optimum length decreases computing runtime. 

Goddard Space Right Center, Greenbelt, Maryland 


One-dimensional Fast Fourier Trans- 
form (FFT) operations work fastest on 
grids whose size is divisible by a power of 
two. Because of this, padding grids (that 
are not already sized to a power of two) 
so that their size is the next highest 
power of two can speed up operations. 
While this works well for one-dimen- 
sional grids, it does not work well for 
two-dimensional grids. 

For a two-dimensional grid, there are 
certain pad sizes that work better than 
others. Therefore, the need exists to 
generalize a strategy for determining 
optimal pad sizes. There are three steps 
in the FFT algorithm. The first is to per- 
form a one-dimensional transform on 
each row in the grid. The second step is 
to transpose the resulting matrix. The 
third step is to perform a one-dimen- 
sional transform on each row in the re- 
sulting grid. Steps one and three both 


benefit from padding the row to the 
next highest power of two, but the sec- 
ond step needs a novel approach. 

An algorithm was developed that struck 
a balance between optimizing the grid 
pad size with prime factors that are small 
(which are optimal for one-dimensional 
operations), and with prime factors that 
are large (which are optimal for two-di- 
mensional operations). This algorithm 
optimizes based on average run times, 
and is not fine-tuned for any specific ap- 
plication. It increases the amount of times 
that processor-requested data is found in 
the set-associative processor cache. Cache 
retrievals are 4-10 times faster than con- 
ventional memory retrievals. 

The tested implementation of the algo- 
rithm resulted in faster execution times 
on all platforms tested, but with varying 
sized grids. This is because various com- 
puter architectures process commands 


differently. The test grid was 512x512. 
Using a 540x540 grid on a Pentium V 
processor, the code ran 30 percent faster. 
On a PowerPC, a 256x256 grid worked 
best. A Core2Duo computer preferred ei- 
ther a 1040x1040 (15 percent faster) or a 
1008x1008 (30 percent faster) grid. 

There are many industries that can 
benefit from this algorithm, including 
optics, image-processing, signal-process- 
ing, and engineering applications. 

This work was done by Bruce H. Dean, 
David L. Aronstein, and Jeffery S. Smith of 
Goddard Space Flight Center. Further infor- 
mation is con tained in a TSP (see page 1 ). 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries 
concerning nonexclusive or exclusive license 
for its commercial development should be ad- 
dressed to the Patent Counsel, Goddard Space 
Flight Center, (301 ) 286-7351. Refer to GSC- 
15678-1. 


© Spatial Query for Planetary Data 

This technology is extensible to Earth science and satellite monitoring and surveillance. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Science investigators need to quickly 
and effectively assess past observations of 
specific locations on a planetary surface. 
This innovation involves a location- 
based search technology that was 
adapted and applied to planetary sci- 
ence data to support a spatial query ca- 
pability for mission operations software. 

Conventional databases of planetary 
datasets are indexed and searchable by 
various metadata, such as acquisition time, 


phase of mission, and target. Searching 
these datasets will produce enormous 
datasets that are difficult, or impractical, 
to browse through to identify observations 
of very specific targets. For queries at spe- 
cific locations, it is fundamentally more ef- 
ficient to specify the location as the target 
of the query; and to have the database 
search based on the location of the data 
rather than metadata that is only indirectly 
or tangentially related to location. 


High-performance location-based 
searching requires the use of spatial data 
structures for database organization. Spa- 
tial data structures are designed to organ- 
ize datasets based on their coordinates in 
a way that is optimized for location-based 
retrieval. The particular spatial data struc- 
ture that was adapted for planetary data 
search is the R+ tree. The R+ tree 
arranges data as a set of nodes that repre- 
sents bounding rectangles. Every leaf 
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